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Abstract. - Slow crack growth in a model of homogenous brittle elastic material is described 
as a thermal activation process where stress fluctuations allow to overcome a breaking threshold 
through a series of irreversible steps. We study the case of a single crack in a flat sheet for which 
analytical predictions can be made, and compare them with results from the equivalent problem 
of a 2D spring network. Good statistical agreement is obtained for the crack growth profile 
and final rupture time. The specific scaling of the energy barrier with stress intensity factor 
appears as a consequence of irreversibility. In addition, the model brings out a characteristic 
growth length whose physical meaning could be tested experimentally. 



Introduction. - Although tensile rupture of atomic bonds requires a stress comparable 
to the Young's modulus, brittle solids commonly rupture at a much lower applied stress 
(typically 3 orders of magnitude lower). Griffith's pioneering work pP has clarified the origin 
of this apparent weakening, postulating that small cracks usually preexist in real solids, with 
stress concentration at the crack tip strongly enhancing rupture. A somewhat similar and 
quite striking effect is the occurrence of failure even though the solid is stressed below its 
experimental breaking threshold (i.e., even if stress concentration due to flaws is taken into 
account). The physical process, sometimes referred to as subcritical rupture, leads to a delay 
in the time for complete macroscopic failure of the solid, with a strong dependence on the 
amplitude of the applied stress. 

A possible driving mechanism for subcritical damaging processes is thermal activation as 
supported by the early experimental work of Brenner and Zhurkov j2JEl ■ Zhurkov introduced 
a kinetic concept of strength of solids, where time to rupture follows an Arrhenius law with an 
energy barrier decreasing with increasing stress fnterestingly, there is still debate about 
whether temperature fluctuations might be sufficient or not to nucleate microcracks and trigger 
crack growth. Recent theoretical works 300121 have emphasized the effect of disorder in 
decreasing the effective energy barrier (conversely increasing the "effective temperature"). 
Other authors have used equilibrium statistical thermodynamics to study how cracks 

might naturally appear from thermal fluctuations in an otherwise homogenous material which 
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goes into a metastable state when stretched. However, another fundamental ingredient not 
usually taken into account is the irreversibility of the rupture mechanism, leading to dynamics 
out of equilibrium. 

Temperature fluctuations in real materials are usually considered too weak to be able 
to overcome a potential energy barrier estimated as the free energy cost to reach Griffith's 
critical crack length. But, the energy cost would be much smaller if rupture was considered 
as an irreversible process, where thermal fluctuations need only to be able to break atomic 
bonds one after another. Irreversibility has already been taken into account by Golubovic et 
al. |10|. with the introduction of a minimum crack opening beyond which the crack can not 
heal back. However, Golubovic introduces a complex behaviour of stepwise growth involving 
several bonds at a time, and eventually find that lifetime is dominated by a critical crack 
length in the same fashion as predicted directly from Griffith's criteria, albeit with a smaller 
critical length and a different power law dependence on applied stress. In this paper, we 
present a different approach by choosing ab initio a distribution of stress fluctuations directly 
linked to thermal fluctuations. Since our goal is to describe irreversible crack growth, we are 
considering the general case of a preexisting initial crack in a flat sheet of material. In order to 
clarify our approach, we derive first the energy barrier corresponding to our geometry in the 
framework of Griffith's energy concept. Then, we give analytical solutions of our model which 
yield the complete growth dynamics for the case of a single crack. Finally, we compare the 
predictions of our model to the results obtained with a numerical simulation of a bidimensional 
elastic system in antiplane deformation. 

Energy barrier from Griffith's energy concept. - Griffith's prediction of a critical crack 
size beyond which there is rupture, i.e. irreversible and fast crack growth, is derived from a 
potential energy taking into account elastic energy due to applied stress a and surface energy 
7 needed to open a crack as a function of a unique order parameter, the crack length £. For a 
bidimensional geometry consisting of a flat sheet with a crack perpendicular to the direction 
of stress, the potential energy per unit thickness of the sheet reads: 



where Y is the Young modulus and Ug is the elastic energy in the absence of crack. This 
energy reaches a maximum for a critical crack length £ c beyond which no stable state exist 
except the separation of the solid in two broken pieces. From this expression of the energy, it 
is clear that a stressed solid without crack is in a metastable state with a lifetime depending 
on the possibility to nucleate a crack with critical length. 

Many authors p HTHHT] have used models essentially inspired by Griffith's energy concept 
to show that the lifetime should follow an Arrhenius law with an energy barrier scaling as 
AU ~ er~ 2 in a bidimensional geometry and AU ~ er~ 4 in a three-dimensional one. Following 
these authors, but taking into account the preexistence of a stable initial crack of length £i 
and introducing stress intensity factors Ki = a^J^li and K c = o-J\£ c , the energy barrier 
becomes : 
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Note that the energy barrier is a linear function of the initial crack length multiplied by a 
function depending only on stress intensity factors. This choice for the energy barrier implicitly 
assumes that there is a possibility for the crack to explore reversible states of crack length 
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between initial and critical crack length, and cannot be a good description if irreversibility of 
crack growth is to be taken into account. 

Energy barrier for an irreversible crack growth. - Let us reflect on the fact that thermal 
fluctuations induce stress fluctuations in the material. The uniaxial loading state of an ho- 
mogeneous solid at fixed temperature is described by its free energy density: <p(o-) = a 2 /2Y. 
Treating stress as a fluctuating internal variable in a fixed volume V , the probability to find 
a given stress is proportional to a Boltzmann factor exp(— ipV/kT). Expanding free energy 
about an equilibrium position cr, the distribution of stress o~f is : 



p((Tf) 
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where k is Boltzmann constant, T temperature and (Act) 2 = kT/(Vd 2 ip/da 2 ) = kTY/V [T2"]. 
When a crack is present, stress concentration increases the probability that breaking occurs 
at the crack tip rather than anywhere else. We assume that stress distribution at the crack 
tip remains the same as eq. © despite the strong divergence of stress and the breakdown of 
linear elasticity. Since the stress intensity factor K w ay/I gives a measure of stress intensity 
close to the crack tip for a crack with length measurement £ when external load is a, we choose 
to work directly with this quantity. The threshold for rupture at the crack tip will be given 
by a critical value of stress intensity factor K c as is usual laboratory practice. The cumulative 
probability that stress intensity is larger than a given value r/ is then : P{rj) = J p (x) dx. 
Since breaking is assumed to be an irreversible process, the typical velocity of crack should 
be directly proportional to the probability to have fluctuations giving stress intensity larger 
than rj = K c : 
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where E is a dimensional constant proportional to the Young modulus. In eq. the last 
equality is valid as long as kT <C rj 2 /2E. We also introduced a typical velocity Vq which 
represents the crack velocity when the condition for crack advance is verified at all times 
(P = I). This quantity is typically the ratio of a microsopic scale (interatomic distance) and 
a characteristic time (inverse vibrational frequency). Since if is a function of crack length £ , 
eq. Q is a differential equation for crack evolution. To solve this equation requires additional 
approximations since the dependence of stress intensity factor on crack length is non-linear. 
First, we introduce a reduced crack length <f> = (£ — £,)/(l c — £i) to measure crack evolution 
as it grows from its initial value (</> = 0) to its ultimate stable value {(f) — I). Then, stress 
intensity factor can be written: 



K ~ o^J~£~ — ay/U + (4 - £i)cf> ~ Ki 



1+ 2U 



(5) 



where the last equality is a reasonable approximation giving less than a 2% error on stress 
intensity factor as long as cj) < 1/2 and £ c < 2£i. Another approximation will be to take K — Ki 
in the pre-factor of the exponential, because neglecting the variation in stress intensity factor 
leads only to a logarithmic correction of the crack velocity. As a consequence of the last 
approximation, the crack velocity will tend to be underestimated. Solution of the differential 
equation (0J is then : 

t = t \l -exp ( -— H (6) 
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where r gives the lifetime of the sample before fast rupture: 



T 



T exp 



(K c - Kjf 
2EkT 



(7) 



with 




(8) 



and <p c is related to a characteristic growth length A: 



A 2EkT li 



(9) 



Note that the crack velocity : 41 = A/(r — t), diverges as time comes closer to lifetime r, 
which simply means that when time r is reached slow crack growth due to thermal activation 
is no longer the driving mechanism, and a crossover towards fast dynamic crack propagation 
will occur. The lifetime r appearing in eq. Q follows an Arrhenius law with an energy barrier 
function of initial and critical stress intensity factors similar to eq. but instead, there is 
no additional proportionality to initial crack length. A similar scaling for the energy barrier 
was found by Marder |13| . 

Results from the simulation of a 2D elastic spring network. - In order to check predic- 
tions from the previous analysis, we model a bidimensional elastic system as a network of 
springs forming a square lattice whose nodes can move only along an axis perpendicular to 
the undeformed plane of springs. The elastic restoring force of the spring is proportional to 
the variation in displacement along the moving axis. This is a simplified model of an antiplane 
deformation. A constant force is applied at two opposite sides of the lattice, the direction of 
the force being parallel to the moving axis, but reversed from one side to the other. Starting 
from an equilibrium configuration, an external force obtained from a normal distribution is 
applied in parallel to each spring. A new static equilibrium configuration is determined to find 
the fluctuations of spring forces. In this procedure, fluctuations are quasi-statically coupled 
to a temperature bath. Making the clastic constant between force and displacement equal to 
unity, the temperature coefficient kT is numerically identical to the variance of spring force 
fluctuations. Whenever the force on a spring exceeds a breaking threshold f c , the spring is 
cut and, since it is never repaired, the process is irreversible. The crack itself is modelled 
as a series of broken parallel springs, with the crack direction parallel to the sides where the 
constant force is applied. The size of the square lattice (100 x 100 springs) is determined to 
reduce finite size effects and obtain the correct scaling of stress intensity factor with applied 
stress and crack length. Time scale in the simulation is given as the constant time between 
two configurations of force fluctuations in the system, and length scale as the distance between 
two nodes of the lattice. 

The distribution of lifetimes r obtained in the simulations is a decreasing exponential with 
a typical width yj (At 2 ) /t w 0.5. In fig. ^ an d HI we plot the lifetimes averaged over an 
ensemble of 10 to 50 simulations for samples with different initial crack lengths and different 
temperatures. In fig. ^ the logarithm of the lifetime is plotted as a function of the energy 
barrier given by eq. J5J) and scaled by temperature coefficient kT. For a given initial crack 
length, a good scaling is obtained at various temperatures. However, when the initial crack 
length is changed, the scaling with the energy barrier does not work at all. In fig. the 
logarithm of the lifetime scaled by To is plotted as a function of the energy barrier given by 
eq. (0. This energy barrier scaling appears to be correct whatever is the initial crack length or 
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temperature. It seems that the lack of dependence on initial crack length in the energy barrier 
given by eq. J7J is a signature of irreversibility. Indeed, if one considers the potential energy 
landscape obtained from eq. I[T|). we can describe a growth process where small irreversible 
steps 8 occur, with the energy barrier to overcome each of these steps proportional to the 
derivative of the potential energy. It can be seen from eq. Q that this energy barrier would 
be expressed only as a function of stress intensity factor as in eq. (jjj. For each step, the 
energy barrier can be written AU = ^jj<5 = (K% — K 2 } S/Y, and the waiting time for the 
crack to advance with step 6 is given by an Arrhcnius law. Integrating the time it takes to go 
from £i to l c gives in the limit 8 <C 1% the following scaling relation: 



Although this approach does correct the problem of scaling with £j, it completely fails to give 
the correct scaling of the energy barrier with stress intensity factors; furthermore, for a given 
initial crack length, temperature scaling is lost as can be seen from the insert in fig. [21 

To determine the characteristic growth length A which appears in our model, it is necessary 
to look at the complete growth dynamic. As for lifetimes, there is a strong dispersion in growth 
curves obtained for different simulations with the same temperature and initial crack length. 
In order to get statistical information about crack growth, we determine the average time it 
takes for the crack to reach a certain length. In fig. El which shows flowing time as a function 
of crack length, the exponential behaviour obtained in eq. © is recovered and yields a value 
of A. In fig. 01 we plot A scaled by the coefficient a — 2E£i/[Ki(K c — Ki)] as a function of kT. 
The linear dependence of A on temperature predicted by eq. © appears clearly for various 
initial crack lengths and applied forces. In addition, the correct rescaling of data with initial 
crack lengths is mainly related to the linear dependence of A on t Dispersion of data for a 
given temperature and initial crack length occurs systematically above prediction of the model 
(see solid curve in fig.@J|. This can be understood as a consequence of the approximations in 
our model underestimating crack velocity. This interpretation is consistent with the increase 
in dispersion of the data as £ c increases. 

Conclusion. We have shown that by taking into account distribution of stress fluc- 
tuations, we can describe thermally activated and irreversible crack growth and obtain the 
correct scaling of lifetime with applied stress and initial crack length. Previous predictions 
based on Griffith's potential energy [§1E31E] fail to take into account irreversibility of the 
rupture process, thus do not apply to our configuration. Our challenge is to find a way to 
include irreversibility in a problem of total potential energy minimization, since simple ar- 
guments fail to give the correct behaviour. More work should be done to check that stress 
divergence at the crack tip does not preclude the derivation of stress distribution presented 
in this paper. Direct comparisons with experiments can be performed to test the relevance of 
the characteristic length scale for crack growth introduced in our model, as well as scaling of 
both lifetimes and characteristic growth length on initial crack length. 




(10) 
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A U/kT 



Fig. 1 - Logarithm of lifetime as a function of the energy barrier as predicted by eq. Failure of 
data to scale with initial crack length is the main observation. Straight lines are a guide for the eye. 




Fig. 2 - Logarithm of lifetime as a function of the energy barrier as predicted by eq. Q. Scaling 
works for various temperature and initial crack length. The straight line, slope 1, is prediction from 
eq. Q. The inset shows failure of scaling with energy barrier as predicted by eq. IIUI . 
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Fig. 3 - Rescaled time as a function of rescaled crack length for various initial crack lengths and 
temperatures. For each set of data, A is extracted as the only fit parameter since lifetime is known. 
Solid line corresponds to eq. ©. 
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Fig. 4 - Linear dependence of A/a with temperature coefficient kT for various initial crack lengths 
and applied forces. The solid line shows the behaviour expected from the model (slope = 1). 



